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Abstract 

The complex fragments (or intermediate mass fragments) observed in the 
low-energy 

58 Ni+ 58 Ni ^ii6 Ba * reac tion, are studied within the dynamical clus- 



ter decay model for s-wave with the use of the temperature-dependent liquid 
drop, Coulomb and proximity energies. The important result is that, due to 
the temperature effects in liquid drop energy, the explicit preference for a-like 
fragments is washed out, though the 12 C (or the complementary 104 Sn) decay is 
still predicted to be one of the most probable a-nucleus decay for this reaction. 
The production rates for non-a like intermediate mass fragments (IMFs) are 
now higher and the light particle production is shown to accompany the IMFs 
at all incident energies, without involving any statistical evaporation process in 
the model. The comparisons between the experimental data and the (s-wave) 
calculations for IMFs production cross sections are rather satisfactory and the 
contributions from other £-waves need to be added for a further improvement of 
these comparisons and for calculations of the total kinetic energies of fragments. 

Pacs No. 25.70. Jj, 23.60.+e, 23.70.+j, 24.10.-i 
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1 INTRODUCTION 



In low-energy heavy-ion reactions, intermediate mass fragments (IMFs) i.e. complex 
fragments with masses lighter than A<20, are emitted from both the light (A~60) 
and medium-mass (A~110) compound systems, and may arise as multiple "clusters", 
and are accompanied by multiple light particles (Z<2) production. This light-particle 
emission, associated with evaporation residues production, is very satisfactorily un- 
derstood as the fusion-evaporation process from an equiliberated compound nucleus 
in the statistical Hauser-Feshbach formalism [1, 2, 3, 4, 5], described in evaporation 
codes such as LILITA and CASCADE or many other codes. The Hauser-Feshbach 
formalism has been extended to include the above noted observed complex IMFs, for 
example in the BUSCO code [2] or in the Extended Hauser-Feshbach Method based 
on the scission-point picture [5] . An alternative process for the IMFs production is the 
binary decay of a compound nucleus in the statistical fission model of Moretto [6] or of 
Vandenbosch and Huizenga [7]. The fission model of [6] is used in the GEMINI code [1] 
for the heavier compound systems (A>100), whereas the saddle-point transition-state 
model based on fission model of [7] is used for the lighter systems (A<80) [4]. The 
emission of light particles (n, p, a or 7-ray) in both of these fission models is still 
calculated within the Hauser-Feshbach formalism [8]. 

We know that the fission and cluster-decay are competing processes with the latter 
being more prevelent for lighter fragments [9]. In view of this result, in a very recent 
paper [10], we have proposed for the first time an alternative new "cluster decay" 
process for the IMFs emitted from an excited 116 Ba* compound nucleus produced in 
the low-energy 58 Ni+ 58 Ni reaction. The IMFs are shown to arise as multiple " clusters" 
of masses A<20, and only at E^ >200 MeV, in agreement with experimental results of 
Beckerman et al. [11]. This result has its origin in the macroscopic liquid drop energy, 
since the shell effects are almost zero at the excitation energies involved. The multiple 
light particle (Z<2) emission, other than via the statistical evaporation process (not 
treated in this model), is also shown to become equally probable at higher energies 
where only the pure liquid drop energies enter the calculations. The macroscopic 
liquid drop energy was, however, taken to be temperature-independent [10] and only 
the angular momentum t = case was considered. Also, a generalized liquid-drop 
model has been proposed recently by Royer and Moustabchir [12] to describe the light- 
nucleus emission. 

In this paper, we extend our previous work [10] by using the temperature-dependent 
liquid drop energy (Vldm(T)) from Davidson et al. [13], since the compound nucleus is 
hot enough at the excitation energies of interest. The other terms of the potential, that 
constitute the scattering potential V(R), are also taken temperature-dependent, based 
on Refs. [13, 14]. In other words, we use a complete temperature-dependent potential 
and apply it to the medium-mass, hot compound system 116 Ba* formed in 58 Ni+ 58 Ni 
reaction at various (low) incident energies [2, 15, 16, 17]. Also, some of the specific 
data [2, 15] chosen here, and discussed in the following paragraph, was not analysed 
in our earlier work [10]. This system has been of much interest from the point of view 
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of the exotic 12 C (or 100 Sn) cluster radioactivity [17, 18, 19, 20, 21], first expected as 
a ground-state decay [18, 19, 20, 21], but now seems to have observable yields only 
for the decay of the corresponding excited compound system [17]. Another important 
point to note is that 116 Ba is mainly a positive Q-value (Q, ou t) system (negative Q, out 
only for some light decay- fragments). Our model is also applied to a negative Q ou t, 
light compound system 56 Ni [22]. Also, the other, both lighter and heavier, isotopes of 
barium, 112 ~ 120 ' 126 Ba, have been of interest for both the ground-state and excited state 
decay studies [18, 23]. 

The data chosen to test our model calculations in the following are the IMFs pro- 
duction cross sections, measured in (i) the 630 MeV 58 Ni+ 58 Ni— > 116 Ba* reaction [2, 15], 
whose kinematical analysis of fragments (the similarity of Gaussian shaped energy spec- 
tra of fragments in coincidence with 7 rays, and the flat angular distribution of their 
ground and excited states) revealed that the IMFs with 4<Z<9 are consistent with 
the decay of a compound nucleus formed in a complete fusion reaction and that the 
effective temperature of the emitting system, estimated on the basis of a simple sta- 
tistical approach for primary fragments, is on the average about 2.2 MeV, which is 
about half of the compound nucleus temperature T (=4.42 MeV, based on Eq. (13) 
below); and (ii) the 348, 371 and 394 MeV 58 Ni+ 58 Ni reaction [17], where the pure 
Hauser-Feshbach BUSCO predictions do not match the measured cross sections (the 
same code is used with a reasonable success for the above mentioned data at 630 MeV 
[2, 15]). In other words, so far only the pure Hauser-Feshbach analysis is made for both 
the experiments. Also, like in [2, 15], the total excitation energy (TXE) in the exit 
channel is observed to be too small (~50 MeV), compared to the available compound 
nucleus excitation energy E^ N (=96.4 to 130.9 MeV). As is already discussed in our 
earlier work [10], such data on TXE (or TKE) are best treated in a fission or a cluster 
decay model, and we use here the cluster decay model. 

The dynamical cluster-decay model for hot compound systems, a reformulation of 
the preformed cluster model (PCM) of Gupta and collaborators [24, 25, 26, 27, 28] for 
ground-state decays, is presented in Sec. 2 and its application to the hot 116 Ba* nucleus 
data is made in Sec. 3. The experimental data for IMFs production cross sections are 
taken from Fig. 6 of [2], and from Table 1 or Fig. 2 of [17]. It may be noted that 
any discussion of the emitted light particles in this paper is of those predicted here to 
occur at higher energies, and for the £=0 case only. The (statistical) evaporation of 
light particles, that may occur just before the beginning of the binary decay process 
of cluster emission studied here, is not included in this paper. Finally, a summary and 
discussion of our results is added in Sec. 4. 



2 THE DYNAMICAL CLUSTER DECAY MODEL 

The dynamical cluster-decay model developed here for hot compound nucleus is an 
adaptation of the preformed cluster model (PCM) of Gupta et al. [24, 25, 26, 27, 28] 
for ground-state decays, which itself is based on the well-known quantum mechanical 
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fragmentation theory, the QMFT [29, 30, 31], given for fission and heavy-ion reactions 
and used later for the prediction of exotic cluster radioactivity [9, 32, 33]. In the 
preformed cluster model, the decay half-life Ti, or the decay constant A, is defined as 

A = ^ = W, (1) 

2 

where Po, the preformation probability, refers to the motion in the mass asymmetry 
coordinate 77=(Ai-A2)/(Ai+A2) and P, the penetrability, to the motion in the relative 
separation coordinate R. This means that, like in PCM, we treat the distribution of 
the IMFs as a dynamical collective mass motion of preformed fragments or clusters 
through the barrier. In principle, the rj and R motions are coupled, but, in view of 
the defining Eq. (1) and our earlier works [29, 30, 34], these are taken to be decoupled 
here. The v$ is the assault frequency, defined below. In terms of the partial waves, the 
compound nucleus decay or the fragments production cross section 



a = £E(2£+i)P P; k = f^ (2) 

with \i = [AiA2/(Ai + A 2 )]m = \Am{\ — rj 2 ) as the reduced mass and £ c , the critical 
(maximum) angular momentum, defined later, m is the nucleon mass. This means that 
A in (1) gives the cross section for £ = case, with a normalization constant z/ , instead 
of the ir/k 2 in (2). However, in the present calculations, made for £ = (s-wave) case, 
the normalization constant is obtained empirically from the experimental data. 
The preformation probability 

P = V^IV^(A))| 2 (2A4), (3) 

(i=l or 2), where ip{rj) are related to ip u (i]), u— 0,1,2,3... (see Eq.(12)), the eigen- 
solutions of stationary Schrodinger equation in 77, at fixed R, 

f -iSI +li, « ?l " , ^ f( * (4) 

which is solved at R=R a , the first turning point in the inter- nuclear potential (see Fig. 
1). For the ground-state (T=0) decay R a =Ci+C 2 =C t , fixed empirically, since this 
value of R (instead of R=Ro, the compound nucleus radius) assimilates to a good extent 
the effects of both the deformations $ of the two fragments and the neck formation 
between them [27]. The Q are the Siissmann central radii Cj=Rj-(b/Rj), with the radii 

/q /q 

Rj=1.28A/ -0.76-l-0.8Aj fm and surface thickness parameter b=0.99 fm. Thus, by 
raising the first turning point R a from R a =R to R a =C t or, in general, —C t +J2 8R(Pi), 
the deformation effects of the two fragments (and the neck formation between them) 
are included here in the scattering potential V(R,T=0), within the extended model of 
Gupta and collaborators for both the light and heavy nuclear systems [27, 35, 36]. The 
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neck-length parameter ^8R{f3i), introduced at the scission configuration, simulates 
the two-centre nuclear shape parametrization and is equivalent to the lowering of the 
barrier for the case of deformed fragments [27] . This method of inclusion of fragment 
deformations in the present calculations is quite similar to what has been achieved in 
both the saddle-point transition-state model of Sanders [3, 4] and the scission-point 
Extended Hauser-Feshbach Method of Matsuse and collaborators [5]. The work of Ref. 
[27] also shows that the alternative of calculating the fragmentation potential V(rj) and 
scattering potential V(R) for actual deformations of the nuclei is not practical since the 
experimental deformation parameters for all the possible fragments (Ai,A 2 ), required 
for calculating V(f]), are generally not available. Note that the deformation effects of 
nuclei in the present calculations are also included via the Siissmann central radii Cj 
and that C t are different for different rj- values and hence C t is C t (r)). The temperature 
effects on radii Rj, and hence on Cj, the neck- length parameter 6R((3i) and the surface 
width b are discussed in the following. 

For the decay of a hot compound nucleus, we follow the prescription of our earlier 
work [10] and define the first turning point 



which means changing R a to R a + AR for the hot compound nucleus. Note that C t 
and AR are also taken temperature-dependent, as defined later. Fig. 1 illustrates 
the scattering potential and other quantities of interest. Apparently, AR depends on 
the total kinetic energy TKE(T), and the corresponding potential V(R a ) acts like an 
effective Q- value, Q e //(T). Since, R a =C t (ri) for T=0, the AR(r)) corresponds to the 
change in TKE at T with respect to its value at T=0 and R a =C t , and hence can be 
estimated exactly for the temperature effects included in the scattering potential V(R). 

We define Q e //(T) in terms of the binding energies of the hot compound system 
and of the cold (T=0) fragments, as 



meaning thereby that the decaying hot compound system, which is at temperature T 
at R=R a , comes out of the barrier and gives in the exit channel two fragments that 
are cold (T=0). These two fragments in the exit channel go to ground state (T=0) by 
emitting some light particle(s) and/or 7-ray(s) of energy (see Fig. 1) 



R a (T) = C t (ri,T) + AR(ri,T), 



(5) 



QeffiT) 



B(T) — [Bi(T = 0) + B 2 (T = 0)] 

TKE(T) 

V(R a ), 



(6) 



E x = B(T)-B(0) 

= Q eff (T)-Q out (T = 0) 
= TKE(T) — TKE(T = 0). 



(7) 



This can be understood by writing from Eq. (7) or (6) 



Q eff (T) = TKE(T) 



5 



= Q out (T = 0) + E x 
= TKE(T = 0) + E. 



(8) 



which means that the fragments are observed in the ground state with Q o «t(T=0) 
(or =TKE(T=0)) and light particle(s) and/or 7-rays emitted with energy E x . The 
remaining excitation energy of the decaying system is then, 

E* CN -E x =\ Q out (T) I +TKE(T = 0) + TXE(T), (9) 

where, in the exit-channel the compound nucleus excitation energy Eq N {—E c _ m _ + Qi n , 
the sum of the entrance-channel centre-of-mass energy E c . m . and the corresponding 
Q-value, Q in =-67.59 MeV for 58 Ni+ 58 Ni^ 116 Ba) plus the (positive) Q-value of exit- 
channel fragments at temperature T, Q, ou t(T), is distributed into the total excitation 
energy (TXE) and total kinetic energy (TKE) of the two outgoing fragments, 

E* CN + QoutiT) = TKE(T) + TXE{T). (10) 

Eq. (9) again shows that the exit channel fragments are obtained with their TKE in 
the ground-state, i.e. with TKE(T=0). The excitation energy TXE(T) in (9) is used 
in the secondary emission of light particles from the primary fragments, which are on 
the average one a-particle to two-nucleons [15], but are not treated here. Instead, 
we compare our calculations with the primary pre-secondary-evaporation fragments 
emission data, if available. 

Note that in Eq. (5), AR depends on 77. In the following, however, based on our 
earlier works [10, 22] and the works of other authors [4, 5] (discussed later), we use a 
constant average value AR, independent of rj, which also takes care of the additional 
J2 SR(Pi) effects of the deformations of fragments and neck formation between them. 
The AR is the only parameter of the model, though we shall see later that the structure 
of the calculated (s-wave) mass spectrum is nearly independent of the exact choice of 
this parameter value. Thus, in our calculations, Eq. (5) is reduced to 

R a (T) = C t (r),T) + AR(T). (11) 

Hence, at temperature T, the preformation factor P in Eq. (3) is calculated for 
R=R a (T) of Eq. (11), with the temperature effects also included in ip(rj) through a 
Boltzmann-like function 

00 

I ^ | 2 = J2 I V I 2 exp(-E u /T). (12) 

The compound nucleus temperature T (in MeV) is given by 

E* CN = (A/9)T 2 -T. (13) 
The penetrability P is the WKB tunnelling probability, 

P = ex V \-\ f Rb {2fi[V(R) - Q ef f}Y' 2 dR}. (14) 

II J R a 

6 



This expression is solved analytically [25]. Here is the second turning point. For 
the penetrability P, Eq.(ll) means that 



V(R a ) = V(C t ( V ,T) + AR(T)) = V(R h ) = Q eff (T) = TKE(T). (15) 

In Eq.(4), the fragmentation potential Vr(?7) at any temperature T is calculated 
within the Strutinsky renormalization procedure, as 

V R ( V ,T) = j2[v LDM (A t ,Z t ,T)]+j2[5U t ]exp(-^) 

i=l i=l I o 

+ 7 ^ + V P (T) + V e (T), (16) 

where, as already stated in the Introduction, the T-dependent liquid drop energy 
Vldm(T) is that of Ref. [13], and shell corrections SU(T) are considered to vanish 
exponentially, with T = 1.5 MeV. The Vldm(T) has the following form, based on the 
semi-empirical mass formula of Seeger [37], 

V LDM (A,Z,T) = a(T)A + P(T)Al + {y(T)-^-) 

// 2 + 2|/K Z 2 / 0.7636 

x 7^ — H r 1 2— 

V A > r (T)A-3 K Zl 

2.29 x , sfrr J(Z,A) 



)+S(T)^p-, (17) 

3 27 Ai 



[r (T)A*]^ A 
with 

I = a a (Z-N), a a =l, 
and, respectively, for even-even, even-odd, and odd-odd nuclei, 

f(Z,A) = (-1,0,1). 

Seeger [37] fitted the constants from the ground-state (T=0) binding energies of some 
488 nuclei available at that time (in 1961), and obtained: 

a(0) = -16.11 MeV, ,3(0) = 20.21 MeV, 

7(0) = 20.65 MeV, 77(0) = 48.00 MeV, 
and from Ref. [38], the pairing energy term 

5(0) = 33.0 MeV. 

These constants apparently need to be re-fitted to the large amount of data on ground- 
state binding energy available now [39], particularly for neutron-rich nuclei. This is 
partly done in Ref. [22] and is extended here to nuclei upto Z=56. The constants 
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required to be fitted again are the bulk constant a(0), working as an overall scaling 
factor, and the proton-neutron asymmetry constant a a , controlling the curvature of the 
experimental parabola, illustrated in Fig. 2 by the open circles for Z=49 nuclei. Note 
the excellent agreement between the present fits (open circles) and the compilation of 
experimental data of Ref. [39] (solid circles). The fits are obtained within 1 to 1.5 MeV 
of the measured binding energies. Table 1 gives as an example the fitted constants for 
1 < Z < 56 nuclei, which includes the ones in Table 1 of Ref. [22]. The kind of 
improvement obtained is further illustrated in Fig. 3 for the fragmentation potentials, 
used in the present calculations (compare the newly fitted binding energies displayed as 
open circles with the experimental binding energies displayed as solid circles [39]). The 
two curves are almost identical, particularly for the positions of maxima and minima. 

The temperature dependence of the constants in Eq. (17) was obtained numerically 
in Ref. [13] by using the available experimental information on excited states of 313 
nuclei in the mass region 22 < A < 250 for determining the partition function Z(A,Z,T) 
of each nucleus in a canonical ensemble (see Ref. [13] for further details) and making a 
least squares fit of the excitation energy E ex (A, Z, T) = Vldm{A, Z, T) — Vldm{A, Z, 0) 
to the ensemble average E ex (A, Z, T) = T 2 ^lnZ(A, Z, T). The a(T), (3(T), 7 (T), rj(T) 
and 5(T) thus obtained are given in Figure 1 of Ref. [13] for T < 4MeV, extrapolated 
linearly for higher temperatures. The 5(T) is constrained to be positive definite at all 
temperatures, with 5(T > 2 MeV) =0. For the bulk constant a(T), the authors also 
give the following empirically fitted equation to a Fermi gas model, 

a(T) = a(0) + ^, 

which is used here in our calculations with the re-fitted a(0). Finally, for the radius 
constant r (T) in (17), the analytical form, taken from Ref. [40], is 

r (T) = 1.07(1 + 0.01T). (18) 

i 

This means that R» (=roA?), and in turn Q, are T-dependent. 

The shell effects 5U, added to the liquid drop energy in Eq. (16), are obtained 
empirically by Myers and Swiatecki [41], for spherical shapes, as 



5U = C 



F(N) + F(Z)_ cA , 
(A/2)* 



(19) 



with 



nx) = ^ i t: 'tot 1 ) {x ~ M *- i} - h x " ~ M ^ ) (20) 




. 5 



where X=N or Z, Mj_i < X < M { . M { are the magic numbers 2,8,14 (or 20), 
28,50,82,126 and 184 for both neutrons and protons. The constants are C=5.8 MeV 
and c=0.26 MeV. We refer to the use of magic numbers 14 or 20 as MS14 or MS20 
parametrization. 
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The Vp in Eq. (16) is the additional attraction due to the nuclear proximity 
potential [42], also taken to be temperature-dependent, as 



V P (R, T) = 47r J R(T) 7 fe(T)$( S , T), (21) 

where -R(T) is the inverse of the root mean square radius of the Gaussian curvature 
and $(s, T) is the universal function, independent of the geometry of the system, given 
by 

<b(«T\ = \ 4(*- 2.54) 2 - 0.0852(s-2.54) 3 for s < 1.2511 
l ' ' \ -3.437e:rp(-^) for s > 1.2511 1 ' 

" a(T) • ( 23 ) 
Furthermore, 7 is the specific nuclear surface tension given by 

'N- Z\ 2 



7 = 0.9517 



1.7826 



MeVfrn' 2 . (24) 



In Eq. (22), s(T) {= R ^p^ ) is the overlap distance, in units of b, between the colliding 
surfaces. The surface width b is [14] 

b{T) = 0.99(1 + 0.009T 2 ). (25) 

Thus, in C, = Ri — (b/Ri), the T-dependence is via both Rj and b. 

The same temperature dependence of r (T), given by Eq. (18), is also used for 
the Coulomb potential E C {T) = Z 1 Z 2 e 2 /r (T)A 1 ^ 3 where the charges Zi are fixed 
by minimizing the potential Vr(i],T) in the charge asymmetry coordinate 7]z={ r L\- 
Z 2 )/(Z 1 +Z 2 ). The Vi is the centrifugal potential. However, like in Refs. [10, 22], we 
consider here only the £=0 case (V^=0) since the s-wave is shown to contain the gross 
information about both the measured fragment emission yields and their total kinetic 
energies. For a detailed description, however, the angular momentum effects must also 
be included, which we wish to do in our next publication. 

The mass parameters B m (rj), representing the kinetic energy part in Eq. (4), are 
the smooth classical hydrodynamical masses [43], since we are dealing here with a 
situation where the shell effects are almost completely washed out. In the classical 
hydrodynamical model, for touching spheres, we use 

B ^ = ^^{ vJX) V (26) 



with 



P{1 } ~ AC t (T){ 2 Ci(T) C 2 (T)) {27) 
v c (T) = nR 2 c (T)C t (T), R C (T) = 0.4C 2 (T), (28) 
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and v t — v± + v 2 , the total conserved volume. Also, C 2 << C\ and R c [^ 0) is the 
radius of a cylinder of length C t , whose existence allows a homogeneous, radial flow 
of mass between the two fragments. Note that here all radial or radial-dependent 
quantities are temperature-dependent, with Ri(T) = r (T)A^ 3 fm in and C t . The 
above expressions are also used for T-independent calculations, but then with Ri = 
1.284 1/3 - 0.76 + 0.84~ 1/3 fm. 

In Eq.(l), the assault frequency u is given simply as 

with the kinetic energy of the lighter fragment E 2 = (Ai/A)Q e ff, where Q e // is shared 
between the two fragments as inverse of their masses. However, for the calculations 
of s-wave cross sections, instead of u , we use an empirically determined normalization 
constant. 

Finally, the temperature-dependent scattering potential V V (R,T), normalized to 
the exit channel binding energy, is 

V V (R,T) = ^^ + V P (T), (30) 

with Ri(T) = r (T)Ay 3 fm. This is illustrated in Fig. 1. For T-independent case 
R { = 1.28A 1 / 3 - 0.76 + 0.8A^ 1/3 fm. Eq. (30) means that for fixed 77, the energies are 
measured w.r.t. -Bi(T) + B 2 (T), and the fragments go to ground state (T — > 0) only 
via the emission of particle(s) or 7-ray(s) of energy E x . 



3 CALCULATIONS 

In this section, we present our calculations in such a way that the role of temperature is 
emphasised in the different terms of the potential. For this purpose, we first look at our 
calculations for the temperature effects only in liquid drop energy and shell corrections, 
i.e. for use of Vldm(T) and 5U(T), but T-independent E c and Vp. Note that this 
is already a first improvement over our earlier calculation in [10], where only shell 
corrections were taken to be T-dependent. Then, as a second and final improvement, 
we include the T-dependence in all terms of the potential, i.e. we use Vldm(T), 5U (T), 
E C (T) and Vp(T). The resulting fragmentation potentials for 116 Ba are given in the 
following two subsections, respectively, for each case. Only s-wave (£ = 0) solutions 
are studied here in the following. 



3.1 Temperature-dependence only in liquid drop energy and 
shell correction 

Figure 4 gives our calculated fragmentation potentials V(^) for cold (T=0) as well as hot 
116 Ba* nucleus at various T-values referring to compound nucleus excitation energies 
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E^^y of the experiments of Refs. [2, 17]. Comparing with our earlier calculation [10], 
here the liquid drop potential Vldm is also T-dependent. The R a - value for each T- value 
is chosen as follows: At T=0, R a =C t , as is already discussed above. For hot compound 
systems, we fix R a =C t + AR with AR = 0.226, 0.252, 0.281, and 0.466 fm (chosen 
arbitrarily), respectively, for T=2.91, 3.06, 3.20, and 4.42 MeV, which correspond to 
the four experimental values of E c m =174, 185.5, 197, and 315 MeV of Refs. [17] and 
[2]. We notice that, as far as the fine structure is concerned, the calculations are nearly 
independent of the exact choice of Ai?-value, which is explicitly illustrated in Fig. 5 for 
T=4.42 MeV. This observation has already been made in our earlier studies on QMFT 
[28, 31]. As a consequence, the same is true for the calculated yields shown below in 
Figs. 8 and 9. Thus, for a given E am ., AR is fixed for all the fragments; on the other 
hand AR was fitted for each fragment separately in our earlier work [10], as well as in 
the statistical model calculations of other authors [4, 5] . Note, the AR in our model is 
the same as the neck length parameter d in Refs. [4, 5]. 

We notice in Fig. 4 that, due to the T-dependence of Vldm, the fragmentation 
potentials show large differences as compared to our earlier calculations in Ref. [10]. 
First of all, the a-like structure, present at all T values in our earlier calculations, 
including the extreme case of T=4.42 MeV where <5U(T) is reduced almost to zero, is 
now washed out, even at the lowest temperature. The potential energy surfaces are now 
smooth for the heavier fragments at all T values, and new non-a fragments appear in 
the lighter mass region. In other words, with increased T, not only the shell structure 
effects go to zero but also the preference for a-like structure, due to the Wigner term, 
vanishes. Also, the light particle (A<3) structure is found to change with T. Thus, with 
T added, the significant structure in V(?7) seems to remain only for a small window of 
light fragments having masses less than ~20, and we refer to it as the "IMF window" 
in the following. 

Figure 6 shows the preformation factor P for the IMF window only, calculated 
for the four chosen experimental temperatures. We notice that the formation yields 
increase as T increases and that some of the fragments are strongly preformed compared 
to their neighbours. Also, the yields for fragments with masses A2 >16 are very small, 
with 12 C being better preformed at higher energies. The lighter particles (A<3), other 
than the statistically evaporated ones (not considered here), are also equally strongly 
preformed and a neutron (n) at lower energies changes to a proton ( 1 H) at higher 
energies. However, the final decay constant (equivalently, the IMF s-wave production 
cross section) is a combined effect of Po and P, the vq being almost constant. The 
result of this calculation is presented in Fig. 7 for all the strongly preformed and 
observed IMFs (Z<9, A<18) at the four experimental energies. Some of the light 
particle emissions (n and 2 H) are not shown here since their Q e ff values are negative 
(or nearly zero) for all the T and R a values considered here (the dependence of the 
decay constant A on R a is already studied in our earlier work [10]; see also Fig. 9 
below). For 1 H, Q e ff is positive only for the two higher energies corresponding to 
T=3.20 and 4.42 MeV. 

The role of the penetrability P is evident in Fig. 7, since some of the strongly 
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preformed fragments, like 4 He and 6 Li, are now shown as less favoured decays (smaller 
decay constants). 4 He, which has the largest preformation factor and decay constant 
for the ground-state decay (see Figs. 3 and 4 in Ref. [10]), lies here lowest for the hot 
116 Ba*. On the other hand, the other two light particle ( 1 H and 3 He) decays become 
the most probable ones at all T values here, and are kind of competing with the so- 
called statistically evaporated particles. Also, the non-a like decays, like 10 B and 14 N, 
with relatively low Po's, are now the most predominant ones. In other words, with the 
T-dependence of Vldm included, the exclusive preference for a-like IMFs seems lost, 
the decay rates for non-a like decays are higher, and the light particle production rates 
are the highest. However, the 12 C decay is still striking since it is the best a-nucleus 
decay even now, particularly at lower energies (inspite of the lower preformation factor 
Po at these energies). At higher energies, the lighter a-like fragment 8 Be competes 
with it, though we know that 8 Be is unbound in the ground state. Note that the 
charge number Z changes with T for some fragments (here for the mass 17 fragment) 
and hence at higher T's the measurements of charge distributions become equally 
important [1, 16]. The sudden increase in the production cross section (A- value) for 
17 F at T>3.2 MeV is due to its enhanced penetrability P, the Q-value effect. (The 
same effect is observed later in Fig. 13 for T-dependence in all terms of the potential, 
but at higher temperature T=4.42 MeV only). 

Fig. 8, upper panel, shows the IMF production cross sections a for Z 2 <9, mea- 
sured as compound nucleus decays at E c . m .=315 MeV [2], compared with our calculated 
s-wave production cross sections i.e. A-values for Ai?=0.466 normalized to the exper- 
imental data for Z 2 =6 fragment decay. The calculated yields are for the energetically 
favoured, most probable fragments in both the mass and charge numbers (see Figs. 4 
and 6) and hence the yields for the same Z 2 but different masses are summed up. For 
example, the fragments contributing to the Z 2 =6 fragment yield are the n C, 12 C, and 
13 C fragments, and the same to the Z 2 =7 fragment yield are the 14 N and 15 N fragments, 
etc. For Z 2 =4, however, we should be cautious in making this comparison between the 
experiment and calculation since 8 Be, being unbound as a ground state, may not be 
measured, but we have included it in our summed yields for Z 2 =4 (the other component 
is 9 Be, with a rather small yield ~ 10~ 4 ). The lower panel shows the same for another 
normalization of our calculations to the Z 2 =7 fragment data. Furthermore, Fig. 9 
presents the same comparison (only for normalization to the Z 2 =6 fragment data) for 
another two AR values. Apparently, the calculations compare reasonably well with 
experiment and are almost independent of the (empirical) normalization constant as 
well as the exact choice of AR value. For the heavier fragments (Z 2 >8), the calculated 
cross sections are rather small firstly because we are considering here the contribution 
of only £ = term and secondly because the temperature effects are included only in 
Vldm and 5XJ. Also, in experiments [2, 15] the contributions from other sources such 
as the peripheral collisions like projectile splitting or deep-inelastic-like mechanism are 
difficult to be separated out for the heavy fragments. Thus, in view of the fact that 
we are dealing here with only the £ = case and that the temperature effects are not 
fully included in the potential, the comparisons between the theory and experiments 
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in Figs. 8 and 9 could be considered rather satisfactory. 

A comparison of our calculations with the experimental data at the other three 
lower energies [17] is shown in Fig. 10 for only the four heavy residues 99 Cd, 100 In, 
101 Sn and 102 In that have been actually measured at each of these three energies, as 
reported in Ref. [17] (see Fig. 2 and Table 1 of [17]). In this case, the calculations have 
been normalized arbitrarily to the Ai=102 fragment. Some calculated fragments show a 
change of fragment charge number since we have considered here only the energetically 
favoured, most probable fragment in both mass and charge numbers which could be 
different from the one actually measured. The comparison shows that the predicted 
yields are in reasonable agreement with the experimental data. 

3.2 Temperature-dependence of the full potential 

Figure 11 gives the calculated fragmentation potentials V(rj, T) for T-dependence added 
in all the terms of the potential (Eq. (16), for i = 0). The T and AR values are the same 
as in Fig. 4 for no T-effects in E c and Vp. (The graph for T=0 is again shown here for 
comparisons). We notice that there is hardly any sizeable effect due to the additional 
T-dependences in E c and Vp, when Fig. 11 is compared with Fig. 4, at least for the 
structure and positions of minima, especially at lower T-values. At higher T-values, 
there are some very small changes, particularly for very light fragments (A 2 <2). It is 
known [44] that Vp does not contribute towards the structure and positions of minima 
in the fragmentation potential and that E c is a smoothly increasing function of 
A 2 . However, the mass parameters B m , though smooth, are now T-dependent, and this 
behaviour seems to enhance the formation probability P , presented in Fig. 12. The 
enhancement is by about one order of magnitude (compare Fig. 12 with Fig. 6 without 
T-dependence in E c , Vp, and B w ). As a consequence, we note that in the formation 
yield of 12 C, relative to its neighbouring fragments, the small peaked structure present 
in Fig. 6 is now absent in Fig. 12. 

Figure 13 presents the result of our calculation for the decay constant A (equiva- 
lently, the s-wave production cross section) as a function of compound nucleus tem- 
perature corresponding to four energies (174, 185.5, 197, and 315 MeV) of the two 
selected experimental data sets of Refs. [2, 17]. We notice important differences with 
respect to the case of no temperature in E c and Vp (compare Fig. 13 with Fig. 7). 
Here, for light particles, instead of only X H and 4 He at higher temperatures, we get 
2 H, 3 H and 4 He at all the temperatures which is favoured by the experimental data 
(for the neutron, the Q e // is still negative). Note that here these light fragments are 
not the statistically treated prompt particles, though the calculated yields are shown 
to be large (calculations made so far are only for the £ = case). Interesting enough, 
12 C decay is still the most favoured a-nucleus decay, though other non-a fragments 
are shown to compete. Also, some fragments having smaller cross sections at lower 
energies are shown to have enhanced cross sections at the higher energy. 

Finally, a comparison of our calculated s-wave cross sections with the measured cross 
sections is made in Figs. 14 and 15. Fig. 14 shows this comparison for the E c . m .=315 
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MeV data [2] with our calculated decay constants (s-wave cross sections) normalized 
to Z 2 =6 fragment data. Similarly, Fig. 15 gives a comparison of our calculations with 
the measured cross sections for the three lower energies of Ref. [17]. Here the data 
are obtained for only the heavy fragments and we have normalized our calculations to 
the mass 102 fragment. The comparisons between the calculations and data are rather 
good at all energies, and better than for the earlier case displayed in Figs. 8, 9 and 
10. The improvement demonstrates the need of including temperature dependence in 
all the terms of the potential. 

4 Summary and discussion 

Summarizing, first the use of T-dependent Vldm is studied for the cluster decay process 
proposed for the IMFs emission [10] induced in heavy-ion collisions at low bombarding 
energies (E c . m . <15 MeV/A). The model is extended and applied to 58 Ni+ 58 Ni^ 116 Ba* 
at various centre-of-mass energies from 174 to 315 MeV. Interesting and important 
enough, with the use of a T-dependent V LDM (for £=0 case only), the explicit preference 
for a-like IMFs (or clusters) is lost, though the 12 C (or 104 Sn) decay (first expected 
from the ground state decay of 116 Ba [18, 19]) is still found to be the most preferred 
a-like fragment. Also, the non-a IMFs and some light particles with Z<2 are shown 
to be produced with larger yields, without the inclusion of any statistical evaporation 
process. With T-dependence included in all the terms of the potential, the production 
yields are enhanced and comparisons with data are improved. Further improvements 
are expected with the inclusion of contributions from all £-waves, which is being worked 
out. 

Also, the TKEs are measurable quantities, but the same are available for the 
58 Ni+ 58 Ni reaction only for E c m .=315 MeV [2] and at only one angle and four charges 
(without the mass identification of IMFs). In our model, as already defined above and 
in Fig. 1, the angle-integrated total kinetic energies are given by TKE=V(R a )=V(Ct + 
AR)=Q e //. Since the calculated production yields are almost independent of the (ar- 
bitrarily) chosen average AR values, only the relative variation and not the absolute 
magnitude of TKEs could be obtained at present in this model. Furthermore, the 
angular momentum contribution (not included as yet) is a must for TKEs calcula- 
tions. Such type of calculations have already been performed [22] for the lighter sys- 
tem 32 S+ 24 Mg— » 56 Ni where the angle integrated TKEs data were available [3]. It will 
be interesting to apply similar considerations to 116 Ba* when similar data on TKEs 
become available. Also, it would be of interest if the statistical fission model, as in- 
troduced in the GEMINI code [1] could also be applied to the 58 Ni+ 58 Ni reaction and 
its predictions compared with the present cluster decay model results as well as to the 
pure Hauser-Feshbach BUSCO results. Furthermore, a similar comparison using the 
alternate picture as generalized liquid-drop model [12] or the so-called unified fission 
models for cluster decay studies [9] is called for since these models do not introduce 
any preformation factor. 
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FIGURE CAPTIONS 

Fig. 1. The s-wave scattering potential, illustrated for 116 Ba* ^ 12 C+ 104 Sn, showing 
the effect of temperature for its inclusion in all terms of the potential (solid line), 
i.e. V(R,T)=E C (T)+V P (T), with the binding energy B(T)=V LDM (T)+5U(T). 
The decay path is shown to begin at R a for the T^O case, along with other 
quantities that are discussed in the text. 

Fig. 2. The mass excess AM (=M j4 -A=NM n +ZM p +B(Z,N)-A in MeV) as a function 
of neutron number N for Z=49 nuclei, calculated by using the 1995 experimental 
data (solid circles) [39], with newly fitted constants (open circles) and with the 
1961 Seeger's constants [37] (open squares). 

Fig. 3. The fragmentation potential for 116 Ba at T=0, using the experimental binding 
energies (solid circles) [39] and the empirical binding energies (open circles) with 
the new constants given in Table 1. 

Fig. 4. The fragmentation potential V(?7,R,T) for the hot 116 Ba* nucleus, calculated at 
the ground-state (T=0, R a =C t ) and at the various nuclear temperatures T, which 
correspond to the bombarding energies of the measurements of Refs. [2, 17], and 
at the indicated R a - values. The T-dependence is included only in Vldm and SU, 
i.e. V(?7,R,T)=Vl j dm(T)-|-(5U(T)-|-E c -|-Vp. The thin vertical lines are drawn to 
mark the exclusive a-like structure present in T=0 case only. 

Fig. 5. Same as for Fig. 4, but for only T=4.42 MeV and at different AR values. 

Fig. 6. The IMFs preformation probability P for the hot 116 Ba* nucleus, calculated 
for the four experimental T-values corresponding to E c . m .=174, 185.5, 197, and 
315 MeV or equivalently T=2.91, 3.06, 3.20 and 4.42 MeV, respectively, using 
the fragmentation potentials displayed in Fig. 4. 

Fig. 7. The cluster decay constant A(s _1 ) as a function of nuclear temperature T for 
emissions of the IMFs and light particles from the hot 116 Ba* nucleus, calculated 
at the four experimental T and R a values used in Fig. 4 that gives the relevant 
fragmentation potentials. The dashed line shows the change of charge number 
of the cluster. As for the experimental data of Refs. [2, 15], only the fragments 
with Z<9 are found to be of a compound nucleus origin. 

Fig. 8. The calculated (s-wave) and measured cross sections cr(mb) for all the IMFs 
produced in decay of 116 Ba* at E c m =315 MeV. The experimental data are taken 
from the Fig. 6 of Ref. [2]. Since in experiments [2] only the charges of fragments 
are measured, in calculations we have summed up the yields for most probable 
fragments with the same charge but different masses. The calculations are made 
for Ai?=0.466 fm and are normalized to experimental data of Z 2 =6 fragment 
(upper panel) and Z 2 =7 fragment (lower panel). 
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Fig. 9. Same as for Fig. 8 (upper panel), but for Ai?=0.578 fm (upper panel) and 
0.8 fm (lower panel). 

Fig. 10. Same as for Fig. 8, but for E c . m =174, 185.5 and 197 MeV [17] and heavy 
fragments, and for normalization with respect to Ai = 102 mass fragment. The 
experimental data are from Table 1 or Fig. 2 of Ref. [17]. Note that the calculated 
and measured Z-values of some fragments are different since in calculations we 
have considered only the energically favoured fragment with most probable charge 
and mass. 

Fig. 11. Same as for Fig. 4, but for the T-dependence included in all terms of the 
potential, i.e. V(r/,R,T)=V LDM (T)+5U(T)+E C (T)+V P (T). Here the potential 
for only the light fragments is shown, since it is symmetrical with respect to 
A 2 =58 (or r]=0). 

Fig. 12. Same as for Fig. 6, but for the potentials given in Fig. 11, i.e. for T- 
dependence included in all terms of the potential. 

Fig. 13. Same as for Fig. 7, but for the potentials given in Fig. 11, i.e. for T- 
dependence included in all terms of the potential. The mass A 2 =l fragment (n, 
in this case) is not shown because Q e ff is negative for the chosen parameters of 
the calculations. 

Fig. 14. Same as for Fig. 8 (upper panel), but for the potentials given in Fig. 11, i.e. 
for T-dependence included in all terms of the potential. 

Fig. 15. Same as for Fig. 10, but for the potentials given in Fig. 11, i.e. for 
T-dependence included in all terms of the potential. 
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Table 1: Re-fitted bulk and asymmetry constants of Seeger's liquid drop energy for 



1< Z < 56. 



z 


N 


a(0) 


a a 


z 


N 


a(0) 


a a 


1 


2 


-15.85 


0.10 


6 


4 


-15.70 


0.10 




3 


-16.95 


0.12 




5,7 


-16.50 


0.10 




4 


-13.00 


0.05 




6 


-16.65 


0.10 




5 


-13.70 


0.12 




8 


-15.90 


0.10 


2 


1 


-15.50 


0.10 




9 


-15.70 


0.10 




2 


-16.00 


0.10 




10 


-15.10 


0.10 




3 


-16.80 


0.30 




11 


-14.80 


0.10 




4,5 


-14.20 


0.30 




12,13,15,16 


-15.00 


0.80 




6 


-13.50 


0.10 




14 


-14.85 


0.80 




7,8 


-13.00 


0.10 


7 


3 


-14.30 


0.20 


3 


1,2,4,5 


-16.60 


0.10 




4 


-15.20 


0.50 




3 


-16.98 


0.98 




5 


-16.20 


0.80 




6 


-13.80 


0.98 




6 


-16.55 


0.80 




7 


-14.30 


0.40 




7 


-16.80 


0.80 




8,9 


-13.20 


0.10 




8 


-16.30 


0.80 


4 


1 


-13.00 


0.01 




9 


-16.20 


0.80 




2 


-14.50 


0.10 




10,11 


-15.90 


0.94 




3 


-16.20 


0.80 




12 


-15.75 


0.94 




4 


-16.98 


0.98 




13 


-15.80 


0.94 




5 


-16.70 


0.60 




14 


-15.65 


0.94 




6 


-15.50 


0.80 




15 


-15.90 


0.94 




7 


-15.30 


0.50 




16 


-16.00 


0.94 




8 


-14.30 


0.10 




17 


-16.10 


0.93 




9 


-14.00 


0.10 


8 


4 


-14.00 


0.94 




10 


-13.30 


0.01 




5 


-15.25 


0.94 


5 


2 


-14.60 


0.10 




6 


-15.90 


0.94 




3 


-16.50 


0.10 




7 


-16.35 


0.94 




4 


-16.60 


0.60 




8 


-16.20 


0.94 




5 


-16.99 


0.10 




9 


-16.18 


0.94 




6 


-16.60 


0.60 




10 


-15.95 


0.94 




7 


-16.30 


0.10 




11 


-15.93 


0.94 




8 


-15.35 


0.10 




12,14 


-15.85 


0.94 




9 


-15.10 


0.10 




13 


-15.90 


0.94 




10 


-14.45 


0.10 




15 


-16.10 


0.94 




11 


-14.10 


0.10 




16 


-16.15 


0.90 




12 


-13.45 


0.10 




17 


-16.30 


0.92 




13 


-13.10 


0.10 




18 


-16.11 


0.92 




14 


-13.00 


0.40 


9 


5 


-15.25 


0.90 


6 


2 


-13.00 


0.10 




6 


-15.90 


0.90 




3 


-13.85 


0.80 




7 


-16.28 


0.90 
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Table 1: Continued 



z 


N 


a(0) 


a a 


Z 


N 


a(0) 




9 


9 


-16.30 


0.90 


21 


15-23,31-38 


-16.42 


0.77 




10 


-16.15 


0.90 




24-30 


-16.38 


0.78 




8,11,17,19,20 


-16.20 


0.90 


22 


16-39 


-16.42 


0.77 




12 


-16.01 


0.90 
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Fig. 1 Balasubramaniam, Kumar, Gupta, Beck and Scheid, 'Emission of intermediate mass fragments from hot 116 Ba 
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Fig. 2 Balasubramaniam, Kumar, Gupta, Beck and Scheid "Emission of intermediate mass fragments from hot 116 Ba " 
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Fig. 3 Balasubramaniam, Kumar, Gupta, Beck and Scheid "Emission of intermediate mass fragments from hot 116 Ba " 
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Fig. 4 Balasubramaniam, Kumar, Gupta, Beck and Scheid "Emission of intermediate mass fragments from hot 116 Ba " 
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Fig. 5 Balasubramaniam, Kumar, Gupta, Beck and Scheid "Emission of intermediate mass fragments from hot 116 Ba. 
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Fig. 6 Balasubramaniam, Kumar, Gupta, Beck and Scheie 1 "Emission of intermediate mass fragments from hot 116 Ba 
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Fig. 7 Balasubramaniam, Kumar, Gupta, Beck and Scheid "Emission of intermediate mass fragments from hot 116 Ba " 
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Fig. 8 Balasubramaniam, Kumar, Gupta, Beck and Scheid "Emission of intermediate mass fragmentsfrom hot 116 Ba " 
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Fig. 9 Balasubramaniam, Kumar, Gupta, Beck and Scheid "Emission of intermediate mass fragments from hot 116 Ba " 
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Fig. 10 Balasubramaniam, Kumar, Gupta, Beck and Scheid "Emission of intermediate mass fragments from hot 11 6Ba. 



100 



75- 



50- 




25- 



- 



. in // 



-25- 



/ \ 



■50- 



1n 




^20Ne 



tiff tl ' ^ J 

V y V 1 22 Na 

T T Wi 
i«E 14N ^ A/ W W ' V 



8Be 



12C 



16Q 



4He 



T 





30 







■ 










■ 


T = 


0.00; 


AR: 


=0.000 


T = 


2.91; 


AR: 


=0.226 


T = 


3.06; 


AR: 


=0.252 


T = 


3.20; 


AR: 


=0.281 


T = 


4.42; 


AR: 


=0.466 


I 

40 


1 I 

50 


1 1 

60 



10 



20 



Light fragment mass number A / 



Fig. 1 1 Balasubramaniam, Kumar, Gupta, Beck and Scheid "Emission of intermediate mass fragments from hot 116 Ba " 
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Fig. 12 Balasubramaniam, Kumar, Gupta, Beck and Scheie 1 "Emission of intermediate mass fragments from hot 11 6Ba 
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Fig. 13 Balasubramaniam, Kumar, Gupta, Beck and Scheid "Emission of intermediate mass fragments from hot 116 Ba " 
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Fig. 14 Balasubramaniam, Kumar, Gupta, Beck and Scheid "Emission of intermediate mass fragments from hot 116 Ba " 
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Fig. 15 Balasubramaniam, Kumar, Gupta, Beck and Scheid "Emission of intermediate mass fragments from hot 116 Ba " 



